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ABSTRACT 

We present a study of galaxy environments to z^2, based on a sample of over 33,000 K-band 
selected galaxies detected in the UKIDSS Ultra Deep Survey (UDS). The combination of in- 
frared depth and area in the UDS allows us to extend previous studies of galaxy environment 
to z > 1 without the strong biases associated with optical galaxy selection. We study the 
environments of galaxies divided by rest frame (U — B) colours, in addition to 'passive' and 
'star-forming' subsets based on template fitting. We find that galaxy colour is strongly corre- 
lated with galaxy overdensity on small scales (< 1 Mpc diameter), with red/passive galaxies 
residing in significantly denser environments than blue/star-forming galaxies to z ~ 1.5. On 
smaller scales (< 0.5 Mpc diameter) we also find a relationship between galaxy luminosity 
and environment, with the most luminous blue galaxies at z ~ 1 inhabiting environments 
comparable to red, passive systems at the same redshift. Monte Carlo simulations demon- 
strate that these conclusions are robust to the uncertainties introduced by photometric redshift 
errors. 

Key words: Galaxy Evolution, Environment, Deep Surveys - Infrared: Galaxies. 
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1 INTRODUCTION 

It has long been known that the properties of galaxies depend on the 
environment in which they are located. Elliptical, non-star-forming 
galaxies occupy more dense regions of space than star-forming, 
disc-dominated galaxies, giving rise to the so -called morphology- 
density relation dOemlerf 19741 ; Dresslerll980l) . The physical origin 
of this relation is still subject to debate, with disagreement mainly 
centering on whether the relation arises due to internal or external 
processes (nature vs nurture). 

Most recent low redshift studies dKauffmann et all 120041 : 
lBaloghetai]|2004h utilise the Sloan Digital Sky Survey (SDSS) 
or the Two-degree-Field Galaxy Redshift Survey (2dFGRS) 
to conduct stat i stical investigations of galaxy environments. 
iKauffmann et al.l J2004t) constrained the specific star formation rate 
(SSFR) using the 4000A break and found that the SSFR (and nu- 
clear activity) depend most strongly on local density, from star- 
forming galaxies at low densities to predominantly inactive systems 
at high densities. 

Studies by Ivan Per Wei et alj d2008h and iBamford et ail 
(2008) found that structure, colour and morphology are mainly de- 
pendent on galaxy mass but that at fixed mass, colour and, to a 
lesser extent, morphology are sensitive to environment. Studies of 
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Ha ( iBalogh et alj|2004h found that it's strength does not depend 
on environment but that the fraction of galaxies with equivalent 
width, Wo(Ha) >4A is environmentally dependant, decreasing 
with increasing density. They also noted that emission line fraction 
appears to depend on both the local environment (~ IMpc) and on 
the large scale structure (~ 5Mpc). 

Studies at higher redshif ts (z ~ 1) h ave used surveys suc h 
as DEEP2 dDavis et alj|2003h and VVDS dLe Fevre et al|l2005h . 



2007) 



DEEP2 investigations i Cooper et alj 120061 Cooper et al 
used the projected third nearest neighbour statistic, studying galaxy 
properties and the colour-density relation respectively. They con- 
cluded that there is a strong dependence on rest frame (U — B) 
colour, with blue galaxies occupying lower density regions but 
showing a strong increase in mean local density with luminosity 
at z ~ 1. This they conclude is consistent with the rapid quenching 
of star formation by AGN or supernova feedback, as ram pressure 
stripping, harassment and tidal interactions, which occur preferen- 
tially in clus t ers, w ould be insufficient to explain these findings. 
ICooper et al] d2007h also observed that the fraction of galaxies on 
the red sequence increases with local density, as in the local Uni- 
verse, but this weakens with redshift and disappears by z ~ 1.3. 

The VIMOS VLT Deep Survey (VVDS ) investigated the red- 
shift and lu minosity evolution of the galaxy colour-de nsity relation 
up to z~1.5 dCucciati et ai]|2 006). In agreement with lCooper et al.l 
(2007) they found that the local colour-density relation progres- 
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sively weakens and possibly reverses in the highest redshift bin 
(1.2<z<1.5). This may imply that quenching of star formation was 
more efficient in high density regions. The VVDS team also ob- 
served that the colour-density relations depend on luminosity and 
found that at fixed luminosity there is a decrease in the number 
of red objects as a function of redshift in high density regions. 
This implied that star formation ends at earlier cosmic epochs for 
mor e luminous/massive galaxies, which is consistent with downsiz- 
ing dCowie et al. 1996). We note, however, that the VVDS survey 
is based on optical /-band selection, and as such will be strongly 
biased against red, passive galaxies at z > 1. Conclusions from 
deep K-selected samples suggest that the galaxy colour bim odal- 
ity is present to at least z~1.5 (e.g. Cirasuolo et alJ 2007b and 
may be still be present at z~2 dCassata et alJ 20081 ; iKriek et al.l 
2008; Willi ams et ai] |2009t) . Furt hermore red galaxies have been 
seen t o strongly cluster at z>1.5 l Daddi et al]l2003l : lOuadri et al.l 
l2007l;lHartlev et alj|2008l; lHartlev et al.ll2010h which suggests that 
a colour-density relation may also exist at these higher redshifts. 

A number of physical processes may be responsible for the 
observed environmental trends. Mergers or tidal interactions can 
tear galactic discs apart and are likely to play an imp ortant role in 
forming the most massive galaxies observed today dToomre et al.l 
ll972l;lFaroukT et alj[l98ll) . Other processes such as gas stripping 
can severely reduce the star formation rate by removing the cold gas 
from galaxies falling into mass ive dark-matter halos dGunn et all 
1 19721 ; Deke l & Birnboim.ll2006f) . Feedback is also t hought to play 
a major role, either from AGN or supernovae (e.g. lBensonll200 5t 



Spri ngeletal.ll2005h . These processes may heat or eject the gas 
within galaxies and thus effectively terminate any further star for- 
mation, which can rapidly lead to the build-up of the galaxy red 
sequence. Finally infalling cold gas in low mass dark matter ha- 
los may fall directly onto the galaxy, whereas in high mass halos 
the gas is thought to be heated by shocks and t herefore remains 
supported dWhite et al.lll99ll ; lBirnboim et alj|2003l) . A key goal of 
observational extragalactic astronomy is to disentangle which of 
these processes are responsible for establishing the bimodal galaxy 
populations observed in the local Universe. 

With the recent advent of deep, wide-field infrared imaging 
in the UKIDSS UDS we can now extend studies of galaxy envi- 
ronments to z>l. Selection in the infrared avoids the major biases 
against dusty and/or evolved stellar populations, allowing us to in- 
vestigate whether correlations observed at low redshift also occur 
at high redshift and how these change over time. The large con- 
tiguous area of this survey also allows us to probe a wide range of 
environments using large samples of galaxies. 

The paper is structured as follows: §2 outlines the data and 
selection criteria used in this work. §3 discusses the method used 
to estimate galaxy environments. The results are then presented in 
§4 and §5, with §6 summarising our conclusions. Throughout this 
paper we assume a ACDM cosmology with f2 m =0.3, f2A=0.7 and 
#0=71 kms" 1 Mpc -1 . 



2 DATA AND SAMPLE SELECTION 

2.1 The UKIDSS Ultra Deep Survey 

This work has been performed using the third data release 
(DR3) of the UKIRT (United Kingdom Infra-Red Telescope) In- 
frared Deep Sky Su rvey, Ultra-Deep Survey (UKIDSS, UDS; 
lLawrence et "alll2007l ; Almaini et al in prep). The UKIDSS project 
consists of 5 sub-surveys of which the UDS is the deepest, with a 
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Figure 1. Rest Frame (U — B) colours for UDS galaxies between 0.75 < 
z < 1.75. Shown as a histogram and as a function of absolute K-band mag- 
nitude. The division between red and blue is shown, using the red-sequence 
slope derived by Hartley et al. (2010). 



target depth of K=25 (AB) ov er a single 4-point ing mosaic of the 
Wide-field camera (WFCAM. ICasali et alj|2007h . giving the UDS 
an area of 0.88 x 0.88 degrees. The 5<r, AB depths within 2" aper- 
tures for the I, H and K-bands are 23.7, 23.5 and 23.7 respectively 
for the DR3, making it the deepest near infrared survey over such a 
large area at the time of release. For details of the stacking proce- 
dure, mosaicing, catalogue extraction and d epth estimation we re - 
fer the reader to Almaini et al. (in prep.) and lFoucaud et al.l d2007t) . 
The field is also covered by deep optical data in the B, V, R, i' 
and z' -bands with depths of Bab=28.4, Vab=27.8, Rab=27.7, 
i4 R =27. 7 and z' An =26.1 from the Subaru-XMM Deep Survey 
(SXDS) dFurusawa et"ai]|2008l 3a, 2" diameter). Data from the 
Spitzer Legacy Program (SpUDS, PFDunlop) reaching 5er depths 
of 24.2 and 24.0 (AB) at 3.6/^m and 4.5/im respectively and U-band 
data from CFHT Megacam (Uab =25. 5; Foucaud et al. in prep) are 
also utilised, which results in a co-incident area of 0.63 deg 2 after 
masking. 

The galaxy sample used here is primarily based on selection in 
the K-band image, upon which we impose a cut at Kab=23 to min- 
imise photometric errors and spurious sources and also to ensure a 
high level of co mpleteness (~100 per cent) and reliable ph otomet- 
ric redshifts (see lCirasuolo e t al. 2010; Hart ley et alj|20101) . Bright 
stars are trivially removed from the combined catalogue by exclud- 
ing objects on a stellar locus defined by 2-arcsec and 3-arcsec aper- 
tures, which is effective to K < 18.1. T he fainter stars are removed 
using a (B - z) - {z - K) diagram dDaddi et alj|2004l) and the 
criterion (z' — K) < 0.3 (B — z) — 0.5. These cuts and care- 
ful masking of bright saturated stars and surrounding contaminated 
regions leave 33,765 galaxies in our sample. 



2.2 Photometric Redshifts 

The magnitudes from the DR3 catalogue were used to determine 
photometric redshifts and stellar ages by x 2 minimisation using 
a large set of templates. This was performed with a code based 
largely on the HYPERZ package (Bolzonella et al. 2000) using 
both average local galaxy SEDs and tem plates from the K20 s urvey. 
Six SEDs of observed starbursts from I Kinney et alJ d_1996) were 
also used to improve the characterisation of young blue galaxies. 
This yeilded photometric redshifts with a 5z/(l+z)=0.008 with a 
standard deviati on of <t=0.034 after the exclusion of outliers (for 
more detail see ICirasuolo et al.l 120071 ; ICirasuolo et alJ l20ld and 
references within). This also provided rest frame magnitudes and 
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colours, of which the rest frame (U—B) and Jf-band absolute mag- 
nitude are utilised here. Sources with an unacceptable fit (x 2 > 15) 
are removed from our sample as these are likely to be unreliable. 
This removes 4% of the galaxy sample, the majority of these are ei- 
ther QSOs (36%), cross-talk (26%) or the minor members of pairs 
or mergers (23%), with the remainder consisting largely of objects 
with very low surface brightness. The fraction of otherwise useful 
objects rejected is therefore <0.6%. 

2.3 Passive Sample 

To define a passive galaxy subset with minimal contamination from 
dusty star forming objects we use a subset of galaxies outlined in 
lHartlev et al.l d2010l) . Templates were used to fit either an instanta- 
neous burst parameterised by an age, or an exponentially decaying 
star-formation rate parameterised by an age and r, the e-folding 
time in the exponentially declining star-formation rate, such that, 

SFR = SFR x e - a9e/T (1) 

where SFR is the star-formation rate at the time of observation 
and SFRo was the initial value. We define a conservative passive 
sample as galaxies that are simultaneously old (age>lGyr) and 
have ongoing star formation with SFR ^ 0.1% of SFRo, and 
a star forming sample with SFR ^ 10% of SFRo, with 3947 and 
22,158 galaxies in each sample respectively. 

To define the red sequence we performed a \ 2 minimisation to 
fit an equation of the form (U — B) = ax Mk + b to the old, burst 
galaxies, definin g the red sa mple to be all galaxies within 3a of this 
fit (see lHartlev et al.l j2010h for a more detailed description). In this 
work, to s eparate red and blue galaxies, we use the red-sequence 
slope from Hartl ey et al.l d2010t) but choose the division between the 
two populations to fit the minimum in the overall colour bimodality 
(as shown in the histogram in Figure[T|l. This leads to a dividing line 
in the colour-magnitude diagram as follows: 

(U - B) = -7.09 x 10" 3 Af x + 0.52 (2) 

This boundary was found to separate the red and blue populations 
effectively to z~1.75. At higher redshift the bimodality in galaxy 
colours is less clear, which may in part be due to photometric er- 
rors. A full examination of this issue and the evolution of the red 
sequence will be presented in Cirasuolo et al. (in prep). Previous 
studies have found evidence fo r an evolution in the l ocation of the 
red sequence with redshift (e.g lBrammer et ai]|2009h . For simplic- 
ity we choose not to model the red sequence in such detail and 
instead use the fixed colour selection boundary given above. We 
note, however, that using an evolving boundary made no signifi- 
cant difference to any of the conclusions presented in this work. 

Table Q] shows the resulting number of red and blue galaxies 
assigned to each photometric redshift bin, including the conserva- 
tive subsamples of passive and actively star-forming galaxies. 



3 ENVIRONMENTAL MEASUREMENT 

We used two methods to calculate galaxy environment: Counts in 
an Aperture and nth Nearest Neighbour. In both methods all the 
galaxies within a photometric redshift bin are collapsed down onto 
a 2D plane and the redshift information within the bin is not utilised 
any further. In the aperture method apertures of IMpc, 500kpc and 
250kpc diameter (physical) are placed on each galaxy and the num- 
ber of galaxies within that aperture are counted (N^ per ). A sample 
of ~100,000 random galaxies are then put down in the unmasked 



regions and the number of randoms within the aperture are counted 
(N^ per ). The number of galaxies within the aperture is then nor- 
malised to give the final density measurement, p/p r - 

jwAper N Tot 

JL = 111 x III (V) 

Pr N, Aper N^ot I ' 

where N^°* and N Tot are then total number of random points and 
galaxies respectively, so that p/p r =l corresponds to a density con- 
sistent with that of a random distribution of galaxies. This method 
was chosen to be the ba sis of this work as it is conceptually sim- 
ple, and as concluded bv lCooper et al.l d2005l) . this technique has a 
distinct advantage in fields masked by a large number of holes. The 
nearest neighbour method would require the exclusion of a large 
fraction of data close to holes and field edges. 

The nth n earest neighbour method was first employed by 
iDresslerl dl98d) . this calculates the distance to the nth nearest 
galaxy, D n in Mpc and is expressed here as a surface density, 

The surface density, £„ is then renormalised such that, 

S n = ^ (5) 

where S is the median density of galaxies within the field. To re- 
duce the effect of the edges, the distance to the nearest edge was 
calculated and if this was less than the distance to the third near- 
est neighbour then the object was removed from the sample. This 
method was only used in this work to test the primary findings of 
the aperture method. The results are presented in the appendix. 



4 RESULTS 

Below we explore the relationship between galaxy colours and en- 
vironment as a function of redshift. Relatively broad redshift bins 
are used to minimise the contamination due to photometric redshift 
errors. These sources of uncertainty are explored further in section 
5. 

4.1 Red and Blue Environments 

Figure|2]shows histograms of galaxy environments within two aper- 
tures of IMpc and 500kpc diameter. These are divided into four 
redshift bins between z = 0.25 and z = 2.25 and into red and blue 
populations. A KS-test was performed on the samples to assess the 
difference between the red and blue populations, with all but the 
final redshift bin showing highly significant differences between 
the two populations. In the first three redshift bins we can exclude 
the null hypothesis that the red and blue galaxies are drawn from 
the same underlying population with a significance in excess of 5a 
significance, falling to <2a and <1<t significance at z > 1.75. 

Figure [3] is a plot of the mean density of red, blue, passive 
(black) and star forming (cyan) galaxies in bins of absolute K-band 
magnitude. Error bars are derived from the error on the mean den- 
sity of galaxies within a given bin. Sources of error are explored 
further in Section 5. The passive and star forming galaxies are de- 
fined in Section 2. As before they are plotted in four redshift bins 
but with an additional 250kpc aperture. This plot illustrates that 
red and/or passive galaxies reside in significantly denser environ- 
ments than blue and/or star-forming galaxies from the present day 
to z ~ 1.5, and this difference is apparent at all luminosities. This 
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Figure 3. The average galaxy overdensity as a function of K-band luminosity, displayed in four redshift bins (top to bottom) and 
using projected apertures of diameter IMpc, 500kpc and 250kpc (three columns). Galaxies are displayed in red, blue, passive and 
star-forming subsets, as defined in Section 2. Note the change in scale for each column. The environments are defined so that p/p r =l 
corresponds to a density consistent with that of a random distribution of galaxies. 





0.25 < z < 0.75 


0.75 < z < 1.25 


1.25 < z < 1.75 


1.75 < 2 < 2.25 


Red 


3,469 


3,500 


2,565 


481 


Blue 


7,367 


7,255 


4,674 


1,711 


Passive 


1,623 


1,395 


722 


87 


Star-forming 


6,938 


6,741 


5,048 


1,817 



Table 1. Table showing the number of galaxies in each sample and in each redshift range, including our conservative subsets of passive 
and actively star-forming galaxies. 



© 201 1 RAS, MNRAS OOO.QJll] 



Galaxy Environments in the UKIDSS Ultra Deep Survey (UDS) 5 




Figure 2. Histograms of the density of galaxies within IMpc and 5()()kpc 
diameter apertures compared to a random sample for red (dashed line) and 
blue (thick line) galaxies. The a values are obtained by performing a KS 
test, representing the significance in rejecting the null-hypothesis that the 
samples are drawn from the same underlying population. 



0.25<z<0.75 1 Mpc 




0.75<z<1.25 1Mpc 



1.25<z<1.75 1Mpc 



1.75<z<2.25 1Mpc 




0.4 0.6 
(U-B) 




i/' 1.5 



0.25<z<0.75 500kpc 




0.75<z<1.25 500kpc 



1.25<z<1.75 500kpc 




-22.5<KMag<-21.5 



-24<KMag<-22.5 



Figure 4. The average galaxy density in a IMpc and 500kpc diameter aper- 
tures compared to a random sample as a function of (U — B) rest frame 
colour. We separate high (black) and low (green) luminosity subsets. Note 
the change in y-axis values between the two different aperture sizes. 



is comp arable to what has been found in the local unive rse by other 
studies dKauffmann et al.l2004 . lvan Per Wei etail2008l) . 

Figure [3] also shows that passive galaxies (shown in black) 
follow a similar density profile to red galaxies but are on aver- 
age in slightly denser environments. This supports the conclusion 
that passive galaxies within the red population are responsible for 
the enhanced environments compared to blue star-forming objects. 
The environments of galaxies that were red but not in the strict 
'passive' sample were also investigated and these were found to lie 
in-between the red and blue galaxy environments, as would be ex- 
pected (these are not shown for clarity). The actively star-forming 
galaxies exhibit the same environmental dependence as the blue 
galaxies, following the same luminosity-density profile. 

In addition to the clear separation of red and blue galaxies, 
we also find a general trend of increasing galaxy density with lu- 
minosity, particularly for blue galaxies and on smaller scales. In- 
specting the two intermediate redshift bins in Figure [3] on scales 
below 500kpc we find that the most luminous blue galaxies appear 
to inhabit environments approaching those of red/pa ssive galax- 
ies. Th ese results are consistent with the findings of ICooper et al.l 
(2007), who observed a strong increase in local density with lumi- 



nosity for blue galaxies at z ~ 1. Our results appear to extend these 
findings to higher redshift, suggesting that the epoch 1 < z < 2 
represents a key period of transformation of massive galaxies from 
the blue cloud onto the passive red sequence. 

4.2 Colour-Density Relation 

In Figure|4]we display average density as a function of (U — B) 
rest-frame colour, with faint (-22. 5< Mk <-21.5) and luminous 
(-24< Mk <-22.5) objects shown in green and black respectively. 
This shows a clear trend for galaxies below z~2, with redder galax- 
ies occupying denser average environments. In the lowest redshift 
bin this change occurs very abruptly at the boundary in the colour 
bimodality at (U — B) reat ~0.7 (see Figure [T). There is also a 
general trend for more luminous galaxies to occupy denser envi- 
ronments. 

In summary, we have found strong evidence that red galax- 
ies are on average found in denser local environments than blue 
galaxies, extending the comparison to a higher redshift than any 
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Figure 5. Showing the density of red and blue galaxies against K-band lu- 
minosity in four quadrants of the UDS field shown in the bottom left, using 
median values of x and y to calculate the position of the centre. The overall 
result for the field at 1.25<z<1.75 is shown in the bottom right with the 
error bars calculated from the error on the mean of the four quadrants, to 
give an estimate of the cosmic variance. 



previous study. ICooper et al] d2007h probed up to z ~ 1.3, finding 
that red and blue galaxies at this redshift occupy indistinguishable 
environments. In contrast we find that there is a significant differ- 
ence between red and blue galaxies between 1.25 < z < 1.75 
(see Figures [2] and [3). This difference is likely to be due to our 
method of selecting galaxies from deep infrared imaging. Very red, 
passive galaxies will have been missed from the R-band selected 
candidates used in the DEEP2 study. This selection effect can be 
quantified using our UDS sample. Considering only the most lu- 
minous galaxies (Ma- ^-22.5) in the redshift range l<z<1.3, if we 
impose the selection criteria of DEEP2 (Ras^24. 1) we find that 
75% of red galaxies would be missed compared to only 20% of the 
blue sample. At fainter luminosities virtually no red/passive galax- 
ies are selected. 



5 INVESTIGATING SOURCES OF ERROR 

In this section we conduct a number of tests to investigate various 
source of systematic uncertainty. We attempt to quantify the cos- 
mic variance in our data by dividing the survey into four quadrants. 
A volume limited study of the brighter galaxies is also performed, 
to allow a clearer comparison between different redshift bins. Fi- 
nally we investigate the effects of photometric redshift errors using 
Monte Carlo simulations. 




-22 -22.5 -23 -23.5 -24 -22 -22.5 -23 -23.5 -24 -22 -22.5 -23 -23.5 
M K (AB) M K (AB) M K (AB) 

Figure 6. As Figure 3, but using only the brightest galaxies with Mj[ ^-22 
(AB) as tracers of galaxy density. Note the change in y-axis values between 
the three different aperture sizes 



5.1 Cosmic Variance 

To attempt to quantify the effects of cosmic variance on our re- 
sults we divide the UDS field into four quadrants. The median val- 
ues of the x and y positions of the galaxies were used to divide 
the field, to ensure that the number of galaxies in each quadrant 
was comparable. The results for the galaxies in the redshift range 
1.25<z<1.75 are shown in Figure [5] The average density versus 
K-band absolute magnitude is shown for the four quadrants, with 
the bottom left panel showing the projected distribution in the field. 
The lower right panel shows the original result from Figure 3 (red 
and blue only) but with the error bars calculated from the error 
on the mean of the densities in the four quadrants in each magni- 
tude bin. It can be seen that the distinction between red and blue 
galaxy environments remains, although there are clearly significant 
variations across the field. We conclude that while field variance 
is clearly present in these data the primary findings remain robust. 
Similar conclusions were drawn from other redshift bins and on 
other scales, which are not shown in the interests of brevity. 

As a cautionary remark, however, we note that strictly speak- 
ing we have only tested the internal field variance, since despite the 
relatively wide field of the UDS (50 x 50 comoving Mpc at z ~ 1) 
we may nevertheless be prone to large-scale cosmic variance due 
to unusual superstructures (Somerville et al. 2004). Testing against 
such effects will require further wide-area infrared surveys, such as 
the forthcoming VIDEO and Ultra- VISTA surveys. 



5.2 Volume-Limited Sample 

Our measurement of environmental density is based on a flux lim- 
ited survey, so by definition we are using fainter galaxies to de- 
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Figure 7. Plot of the reliable spectroscopic redshifts against photomet- 
ric redshifts, highlighting the three highest redshift bins, 0.75<z<1.25, 
1.25<z<1.75 and 1.75<z<2.25, with AGN and radio galaxies removed. 
The red squares and blue circles indicate the red and blue galaxies respec- 
tively, as defined in section 2.2. 



fine environments at low redshift. To investigate the influence of 
this systematic effect we reduced the sample to only the bright- 
est galaxies with Mk ^-22 (AB) and repeated the environmental 
analysis using this volume-limited sample. This allows the galax- 
ies to be used as tracers in the measure of environment over the 
same range in luminosity at all redshifts, to allow a fairer compar- 
ison between epochs. The results are shown in Figure[6] where the 
average density versus K-band luminosity is shown for comparison 
with the flux-limited study in Figure [3] As expected, the volume- 
limited study is much noisier (particularly at low redshift) but the 
results are consistent with the primary findings of Section 4. Red 
and/or passive galaxies are found to occupy denser environments 
on average to z ~ 1.5. 



5.3 Monte Carlo Simulations 

Monte Carlo simulations were conducted to investigate the effect 
of photometric redshift errors on our results, explicitly allowing for 
the effects of outliers and catastrophic errors. Such errors would 
generally dilute differences in galaxy environments, but could po- 
tentially introduce fake overdensities if a large fraction of low- 
redshift structure is incorrectly assigned to a higher redshift bin. 
We note, of course, that the effects of photo-z errors are already 
present in the data, so these simulations can only provide an indi- 
cation of the magnitude of the shifts due to these effects. 

The distribution of spectroscopic versus photometric red- 
shifts used in this field are shown in Figure [TJ with spectroscopic 
redshifts derived fro m a variety of sources dYamada et ai] |2005l : 
Simps on" et alj|2006l ; Akiyama et al., in prep; Simpson et al., in 
prep; Smail et al., in prep) as outlined in Cirasuolo et al. (2010). 
We use the distribution in each redshift bin (higher redshift bins in- 
dicated in Figure UJ to estimate the fraction of galaxies incorrectly 
reassigned from one redshift bin to another. For simplicity two cat- 
egories were adopted: 'catastrophic' errors were defined as those 
with 8z > 0.25, while others are classified as 'normal'. In the ab- 
sence of catastrophic errors we expect > 99.9% of galaxies to show 




Figure 8. Monte Carlo simulation showing the effect of the photo-z errors 
on our measurement of galaxy overdensity, based on 100 simulations. The 
thick black line shows the average of the simulations and the dashed black 
line shows the original result for 0.75<z<1.25 from Figure[5] 
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Figure 9. As Figure[8]but for the 1.25<z<1.75 bin. 



8z < 0.25 at z ~ 1 (Cirasuolo et al. 2010). The resulting fractions 
were then randomly applied to the full sample, shifting galaxies 
into new bins to simulate the effect of photometric redshift errors 
('normal' and 'catastrophic'). 

This was repeated 100 times per redshift bin. Results are 
displayed for the redshift bins 0.75 < z < 1.25 (Fig[8j and 
1.25 < z < 1.75 (Fig(9](. These simulations demonstrate that pho- 
tometric redshift errors tend to dilute the distinction between the 
red and blue populations rather than introduce fake differences. The 
resulting additional sources of error in galaxy density are compara- 
ble to the errors on the mean displayed previously. Only in highest 
redshift bin (z > 1.75) do the simulated populations introduce sig- 
nificant extra scatter (not shown). Given these findings, and the lack 
of spectroscopic redshifts at z ~ 2, we urge caution in interpreting 
our preliminary results in this highest redshift bin. At lower red- 
shifts, however, the environmental differences appear robust. 

Noting that we have only 97 spectroscopic redshifts in the red- 
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Figure 10. The results of Monte Carlo simulations based on the extreme 
scenario in which all incorrectly assigned photometric redshifts lead to blue 
galaxies being incorrectly classified as red, and vice versa. The original 
measures of environment for red galaxies and blue galaxies as a function of 
luminosity are shown as thick red and blue lines. Black dashed (dot-dashed) 
lines show the results of Monte Carlo simulations for the red (blue) popu- 
lations respectively. The corresponding green and magenta lines illustrate 
the influence of doubling and tripling the number of incorrectly assigned 
galaxies. 



shift bin 1.25< z <1.75 (with 64 between 1.25< z photo <1.75), 
where environmental differences appear to be very significant, we 
ran additional simulations using only half of the spectroscopic sam- 
ple. The results were very similar to Figure(9] indicating that we are 
not suffering from low number statistics. 

As an additional test, we conduct further Monte Carlo simula- 
tions to investigate the effects of contamination of the blue galaxies 
by the red population and vice versa. These simulations were mo- 
tivated by the worry that incorrect redshifts may inevitably lead 
to errors in the resulting restframe (U — B) colours. We adopt 
an extreme approach by simulating 100% misclassification, using 
the contaminating fractions outlined above but moving only blue 
galaxies into the red sample and vice versa. We find that these sim- 
ulations have the effect of diluting the environment for both red and 
blue simulations, but once again are insufficient to influence any of 
the results previously presented. To investigate further, the contam- 
inating fractions are artificially increased to double and triple the 
maximal contamination indicated by spectroscopic redshifts. These 
only serve to further reduce the average galaxy densities. Results 
for the 0.75 < z < 1.25 bin are illustrated in Figure [Tol and simi- 
lar results are found for all other redshift bins. 

We conclude that photometric redshift errors (including catas- 
trophic errors) have only served to dilute the differences in envi- 
ronment, but an unrealistic proportion of misclassified galaxies are 
required to significantly influence our primary findings. 



on rest frame (U — B) colours, as well as using 'passive' and 'star- 
forming' galaxies from template fitting. We attempt to quantify the 
effects of cosmic variance, photometric redshift errors and flux- 
limit biases on the resulting environmental measurements. We find 
the following principle results: 

(i) We find a strong relationship between rest-frame (U — B) 
colour and galaxy environment to 2 ~ 1.5, with red galaxies re- 
siding in significantly denser environments than blue galaxies on 
scales below IMpc. These results are robust to the effects of field 
variance, flux-limit biases and photometric redshift errors. The en- 
vironments appear indistinguishable by z ~ 2, but the current lack 
of spectroscopic redshifts at this epoch does not allow a robust test 
of this tantalising signal at present. 

(ii) Selecting 'passive' and 'star-forming' galaxies using tem- 
plate fitting yields consistent results, with the passive subset occu- 
pying slightly denser environments than the global red population 
at all epochs. We conclude that the passive subset of the red galax- 
ies are responsible for the enhanced environments compared to blue 
galaxies. 

(iii) On small scales (< 0.5 Mpc) we find evidence for a posi- 
tive correlation between galaxy K-band luminosity (a good proxy 
for stellar mass) and local density. This trend is particularly clear 
for star-forming and blue galaxies, with the most luminous blue 
galaxies at z ~ 1 — 2 showing average environments comparable to 
passive red galaxies. These findin gs app ear in very good agreement 
with the findings of lCooper et alj (2007) at 2 ~ 1. These results are 
consistent with the identification of high-z galaxies in transition to 
the red sequence in the densest environments. 

To improve on this study, we note that the UDS imaging pro- 
gramme will continue at UKIRT until 2012 and ultimately probe 
~ 1 magnitude deeper than the data used in this work. This will 
be sufficient to detect L* galaxies to 2 ~ 5 in addition to tens 
of thousands of sub-L* galaxies at lower redshift. Furthermore, an 
ongoing ESO Large Programme (UDSz) will soon transform the 
UDS project with the addition of ~ 3000 galaxy spectra. This will 
increase the number of UDS spectroscopic redshifts at 2 > 1.5 by 
an order of magnitude, which we expect to dramatically improve 
the reliability of photometric redshifts and allow us to extend our 
study of galaxy environments and large-scale structure to the cru- 
cial epoch when the galaxy red sequence is first established. 
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6 CONCLUSIONS 

In this paper we have presented a study of the environment of over 
33,000 K-band selected galaxies between 0.25<z<2.25. We have 
used a simple measurement of projected galaxy overdensities to 
study the environments of galaxies divided into red and blue based 
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Figure 11. The nearest neighbour density of red and blue galaxies with the 
nearest neighbour chosen so that the mean distance probed was a diameter 
of IMpc and 500kpc in each redshift bin. 
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APPENDIX: NEAREST NEIGHBOUR METHOD 

In this appendix we use the nth nearest neighbour method as an 
alternative environmental estimator, to provide a test of the re- 
sults m^ejUed_earVierand to allow comparison with previous work 
(e.g. lCooper et al.l2007T) . The nth nearest neighbour was calculated 
for all the galaxies as a function of absolute magnitude, using the 
method outlined in Section 3. The value of n is chosen to probe the 
same mean scale as the apertures method. The results are shown in 
Figure [Til, which appear globally consistent with those derived by 
the aperture method. 
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